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COMPUTATIONAL STUDY OF NEAR-LIMIT PROPAGATION OF 
DETONATION IN HYDROGEN-AIR MIXTURES 


S. Yungster and K. Radhakrishnan 
Institute for Computational Mechanics in Propulsion 
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Abstract 

A computational investigation of the near-limit propaga- 
tion of detonation in lean and rich hydrogen-air mixtures is 
presented. The calculations were carried out over an equiva- 
lence ratio range of 0.4-5. 0, pressures ranging from 0.2 bar 
to 1.0 bar, and ambient initial temperature. The computations 
involved solution of the one-dimensional Euler equations 
with detailed finite -rate chemistry. The numerical method is 
based on a second-order spatially accurate total-varia- 
tion-diminishing (TVD) scheme, and a point implicit, 
first-order-accurate, time marching algorithm. The hydro- 
gen-air combustion was modeled with a 9-species, 19-step 
reaction mechanism. A multi-level, dynamically adaptive 
grid was utilized in order to resolve the structure of the deto- 
nation. The results of the computations indicate that when 
hydrogen concentrations are reduced below certain levels, 
the detonation wave switches from a high-frequency, low 
amplitude oscillation mode to a low frequency mode exhibit- 
ing large fluctuations in the detonation wave speed; that is, a 
“galloping” propagation mode is established. 

Introduction 

Recent experimental and computational results indicate 
that the specific impulse obtained in a hydrogen fueled, 
air-breathing pulsed detonation engine (PDE) is highest 
when it is operated under fuel lean conditions 1 ' 3 . However, 
when the fuel-air equivalence ratio, ([), reached values below 
approximately 0.6, a sharp drop off in thrust and specific 
impulse was observed experimentally, possibly due to hav- 
ing exceeded the detonability limits for hydrogen-air mix- 
tures. 

At fuel-rich conditions, the performance of an air-breath- 
ing PDE is significantly degraded. Nevertheless, a potential 
application of a fuel-rich pulsed detonation rocket engine 
(PDRE) has been proposed in combination with an ejector. 
In such a PDRE-ejector engine, the excess fuel in the PDRE 
tube is mixed and burned with entrained air and subse- 
quently expanded through a nozzle in an unsteady version of 
the ejector-ramjet concept 4 . Unsteady ejectors are poten- 
tially more efficient than their steady-state counterparts 5 . 

There is therefore considerable practical interest in inves- 
tigating the limits of detonability for both rich and lean mix- 
tures. Published work on detonability limits is not extensive. 

Gordon et al . 6 carried out experiments in hydrogen-oxygen 


mixtures with different diluents. Of particular interest to the 
present study are their results for hydrogen-air mixtures 
which exhibited large oscillations in detonation speed near 

the detonability limits. Dupre et al 7 also showed that in 
hydrogen-air mixtures large detonation velocity fluctuations 
were obtained when hydrogen concentrations were lowered 
below a certain level. Because of the large variation in prop- 
agation speed, these marginal detonations are sometimes 
called “galloping” detonations. 

Belles 8 suggested that the detonability limits in hydro- 
gen-air and hydrogen-oxygen systems are the mixture com- 
positions for which the conditions at the von Neumann spike 
lie outside the isothermal branched-chain explosion limits. 
Belles’ analysis included only four chemical reactions. Dove 
and Tribbeck 9 extended the work of Belles by considering a 
more complete combustion mechanism. They showed that 
secondary reactions involving H0 2 cannot be ignored, and 
that the effects of exothermicity are also important and can 
lead to a “chain-thermal” explosion. That is, explosions can 
occur for compositions significantly outside the isothermal 
limit. However, the numerical methods needed for solving 
the stiff ODE’s describing the combustion process were not 
sufficiently developed at the time, and they could not deter- 
mine the boundaries outside the isothermal limits for which 
explosive behavior would still exist. 

In this paper, we investigate first the combustion of 
hydrogen-air mixtures at the von Neumann spike conditions 
using the LSENS code 1011 . Subsequently, we examine in 
more detail the structure and stability of the detonation wave 
for near-limit mixtures by conducting one-dimensional com- 
putational fluid dynamic (CFD) calculations. 

Reaction Kinetics at the von Neumann Spike 

The reaction rate equations were integrated under the 
assumption of a Zeldovich, von Neumann, Doring (ZND) 
detonation. The actual calculation consisted of two parts: 

(i) The determination of the state at the von Neumann 
spike, which is used as the initial condition for the chemical 
kinetic calculation. This state was found using the CEA 
equilibrium code 12 . 

(ii) The integration of the chemical kinetic rate equations 
with the initial conditions found in step one. These calcula- 
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tions were carried out using the LSENS chemical kinetics 
code 11 . 

The combustion mechanism used in the calculations in 
part (ii) and in the CFD computations described later, is 
based on Jachimowski’s model 13 , and is listed in Table 1. It 
consists of 19 elementary reversible reactions among 8 react- 
ing species, with N 2 treated as an inert (i.e., nonreacting) 
species. 

Results are presented in figs. 1 and 2 in the form of igni- 
tion delay times as a function of mixture equivalence ratio 
for two different initial fill pressures: p 0 = 0.4 bar and 
Po = 1.0 bar. Here, the ignition delay is arbitrarily defined as 
the time required to increase the temperature of the gas by 
25 K. For p 0 = 0.4 bar (fig. 1), ignition occurs within a few 
microseconds for the equivalence ratio range 0.6 < <p < 3.0- 
Outside this range, the ignition delay increases very rapidly. 
For p 0 = 1.0 bar (figure 2), a similar trend is observed, but 
the curves are significantly steeper on the lean and rich sides 
of stoichiometric than those at the lower pressure. Also, the 
ignition delay times near stoichiometric conditions are 
smaller at the higher pressure. At <\> = 1, for example, the 
ignition delay is 0.83 jus for pg = 0.4 bar, and 0.36 jus for 
p 0 = 1.0 bar. Figure 2 also includes the value of the deton- 
ability limits computed by Belles 8 as well as experimental 
values obtained by Laffitte 14 , as reported in ref. 8. 

It is clear from figs. 1 and 2 that the limits of detonability 
are probably determined by the steep increase in the ignition 
delay time for very lean and rich mixtures. That is, for mix- 
tures significantly away from stoichiometric conditions, the 
large ignition delay times may prevent the coupling between 
the reaction front and the shock wave needed to form a stable 
detonation wave. Based on this argument, it appears from 
figs. 1 and 2 that the limits of detonability become further 
apart as the pressure is lowered, a result also predicted by 
Belles 8 . This somewhat counterintuitive behavior can be 
explained from the variation of the ignition delay (or explo- 
sion limit) with pressure for hydrogen-air mixtures. Table 2 
gives ignition delay times for a stoichiometric hydrogen-air 
mixture at an initial temperature of 1250 K. It is observed 
that the ignition delay time is not a monotonic function of 
pressure. For pressures between 10 and 40 bars, the ignition 
delay actually increases with pressure. 

The precise value of the maximum ignition delay that 
would allow the establishment of a stable detonation, or the 
changes in the structure and stability of a detonation wave at 
near-limit conditions cannot be determined from these type 
of calculations. To answer these questions a more detailed 
analysis, based on CFD computations that include the 
dynamics of the flow and combustion process (and their 
interaction) is needed. Such analysis is presented in the fol- 
lowing sections. 


CFD Simulations 

As a first step in understanding the behavior of detona- 
tion waves at near-limit conditions, we carried out 
one-dimensional CFD computations. It is well known that 
detonation waves have a three-dimensional structure, and 
therefore, a one-dimensional analysis will not be able to 
investigate the effect of transverse waves. However, such a 
study will be able to examine, in a relatively quick and inex- 
pensive way (particularly when a detailed reaction model is 
used to describe the chemistry) the longitudinal stability 
characteristics of detonation waves. Moreover, a one-dimen- 
sional analysis can serve as a basis for subsequent 
multi-dimensional calculations. 

Since detonation limits are not dependent on transport 

phenomena 6 , we neglect all dissipation effects and therefore 
based our study on the Euler equations. We use the conserva- 
tion form of the unsteady Euler equations and replace the 
global continuity equation with the 9 species conservation 
equations. A detailed description of the set of differential 
equations and additional state and constitutive equations 
needed for system closure are given in ref. 15. 

The numerical method used for solving the governing 
equation set is based on Yee’s second-order spatially accu- 
rate total -variation -diminishing (TVD) scheme 16 , and a 
point implicit, first-order time-accurate marching algorithm. 
This method is a subset of a more general class of BDF 

methods considered by Yungster and Radhakrishnan 17 . It is 
obtained by setting the variables y=0 and (3=1 in equation 4 
of ref. 17, and treating only the chemical source term implic- 
itly. Further details about the numerical algorithm are given 
in ref. 17. 

In order to maintain adequate numerical resolution of the 
detonation wave front, without the need to use hundreds of 
thousands of grid points, a multi-level, dynamically adaptive 
grid is utilized. Figure 3 shows a section of the grid at three 
different times as the detonation wave moves from left to 
right. The grid constantly adapts to keep the detonation front 
within the finest grid level. An arbitrary number of levels can 
be specified. Nine or ten grid levels were used in the present 
study, and 100 points were included in the finest grid level. 


CFD Results 

There are two ways in which a detonation can be formed 
in general: (1) by direct initiation, in which a strong shock 
wave is generated in the tube (with a charge of solid explo- 
sive or by using a high pressure reservoir) and (2) by transi- 
tion from a deflagration. In this paper, we consider only 
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development of detonations with direct initiation. Figure 4 
shows a schematic of the initiation process used in the 
present study. A high pressure, high temperature driver gas, 
consisting of nitrogen, was used in a small region next to the 
head-end of the tube. As the computation is started, a shock 
wave travels to the right and an expansion wave propagates 
to the left, towards the head-end. The shock wave is strong 
enough to initiate chemical reactions in the combustible mix- 
ture. Depending on the mixture properties, the shock wave 
and the combustion front can subsequently merge and form a 
detonation wave. The value of the driver pressure, will 
determine the degree of overdrive of the resulting detonation 
wave. The degree of overdrive /is defined as 

/ = d 2 /d 2 cj 

where D is the actual detonation propagation speed, and D CJ 
is the theoretical Chapman- Jouguet detonation speed. 

The first case considers a stoichiometric mixture of 
hydrogen-air at p 0 = 0.4 bar and T 0 = 298 K and a driver 
pressure ratio, r p = p ( iri\/Po^ of 100. A grid refinement study 
for this case is presented in figure 5. Numerical results are 
presented for four successively refined grids having a mini- 
mum spacing, Ax^, indicated in each figure. The plots 
show the variation in the detonation speed with time. Ini- 
tially the shock wave travels at a constant speed (determined 
by the initial conditions) slightly above 1500 m/s. Chemical 
reactions behind the shock front start generating strong com- 
pression waves (i.e., a reaction shock) after approximately 
10 ps . The reaction shock overtakes the shock front at 
around 13 jus, strengthening it and sharply increasing its 
propagation speed (figs. 5b-5d). Subsequently, the detona- 
tion wave speed decreases gradually. 

It can be observed that the structure of the detonation 
wave changes with grid spacing. On the coarsest grid 
(fig. 5a), the detonation reaches a steady propagation speed 
of D = 2020 m/s after the initial transient phase. For the next 
finer grid (fig. 5b) the computed low frequency oscillation 
appears to be a numerical artifact. With the finest two grids 
used in this study (figs 5c and 5d), the detonation reaches a 
high-frequency, low-amplitude propagation mode; however, 
with the third grid (fig. 5c) a low frequency mode can still be 
observed (in addition to the high-frequency one). The solu- 
tion in Fig. 5d has an average propagation speed of D = 2025 
m/s, nearly identical to the steady solution in fig. 5a. The fre- 
quency of the oscillations, co, in fig. 5d is 0.90 ± 0.05 MHz. 

Figure 6 shows pressure distributions at different times 
for the grid spacing corresponding to fig. 5d. At t = 10.8 ps, 
chemical reactions behind the shock front have produced a 
sharp peak in pressure. No change in the shock front propa- 
gation speeds occurs at this time. At / = 13.7 ps, the reaction 
shock has overtaken the shock front resulting in a spike in 


pressure and propagation speed. After approximately 
t = 38.3 ps, high frequency pressure oscillations become evi- 
dent downstream of the detonation front. 

The theoretical Chapman-Jouguet detonation speed for 
this mixture, as computed from the chemical equilibrium 
code CEA 12 is Dqj- 1956.2 m/s. Therefore, the results pre- 
sented in figs. 5 and 6 are actually those of an overdriven det- 
onation, with a degree of overdrive /= 1.07. The frequency 
of oscillation computed in fig. 5d is in close agreement with 
the results of Sussman 18 , who obtained a frequency of 
1 MHz on a similar piston-driven detonation computation for 
a degree of overdrive of/= 1.1, slightly higher than that con- 
sidered here. 

In the ballistic range experiments of Lehr 19 , which con- 
sisted of spherical projectiles being fired into hydrogen-air 
mixtures, a frequency of oscillation of 1.04 MHz was 
obtained for a projectile speed of 2029 m/s, and for the same 
mixture considered in Fig. 5, but at a slightly higher pressure 
(0.42 bar). 

The high-frequency mode of oscillation can be explained 
by the McVey-Toong wave interaction mechanism 20 . Two 
fundamental processes form the basis of this mechanism: 
(i) when a new reaction front is created, compression waves 
(reaction shocks) are generated that travel upstream and 
downstream from the new reaction front; (ii) when an old 
reaction front is extinguished, it must be accompanied by the 
generation of upstream and downstream rarefaction waves 
which have a strength comparable to the reaction shocks. 
The interaction between these waves and the shock front 
produces the high frequency oscillations, as described in 
detail in refs. 20 and 17. 

To verify that the high-frequency oscillations computed 
in this case are described by this mechanism, the history of 
pressure and density in the region immediately behind the 
detonation front have been plotted in an x-t diagram in fig. 7. 
This figure shows three pulses corresponding to three oscil- 
lation cycles. (The “staircase” appearance of this plot is due 
to the periodic sliding of the grid to keep track of the detona- 
tion front). Careful examination of this figure reveals all of 
the elements of the McVey-Toong mechanism. 

Unless otherwise stated, all subsequent calculations pre- 
sented below, were carried out on the same grid used in fig. 
5d. Figure 8 shows plots of detonation speed as a function of 
time for several lean mixtures. For (() = 0.8, the structure of 
the detonation is very similar to that observed in the stoichio- 
metric case, except that the frequency has decreased to 
0.76 MHz. For 0 = 0.6, the frequency decreases further to 
0.52 MHz and an additional low frequency mode appears. At 
(() = 0.5 the low and high frequency modes are present simul- 
taneously, and at (() = 0.4 the detonation wave almost com- 
pletely switches from a high-frequency, low-amplitude 
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oscillation mode to a low-frequency, high-amplitude one; 
that is, a “galloping” propagation mode is established. 

While there is no generally accepted criterion for defin- 
ing limit behavior, galloping detonations represent the lowest 
propagation mode of detonation in a tube, and it has been 
suggested that the detonability limit should be based on the 
onset of the galloping mode 7 . 

Figure 9 presents the maximum pressure behind the 
shock front as a function of time for the same hydrogen-air 
mixtures as above. This figure shows the increase in the 
amplitude of the oscillations and the switch to the low fre- 
quency mode as the fuel-air equivalence ratio is lowered. 

Figure 10 shows the variation in detonation wave speed 
with time for rich mixtures. Note that as (J) increases, the fre- 
quency of oscillation decreases and the time it takes for the 
detonation to develop also increases significantly. For 
c|) = 4.0, only a few oscillation cycles were captured before 
the detonation wave reached the end of the computational 
domain. The frequencies of oscillation are: 0.62 MHz for 
4> = 2.0, 0.41 MHz for 4> = 3.0, and 0.29 MHz for = 4.0. 

Pressure effect 

The effect of pressure on detonation wave structure was 
examined and the results are presented in figs. 1 1 and 12, for 
p 0 =0.2 and 1 .0 bar respectively. 

The theoretical Chapman- Jouguet detonation speed for a 
stoichiometric hydrogen-air mixture at p 0 = 0.2 bar and 
7*0 = 298 K is 1939.5 m/s. The average detonation speed 
computed in fig. 1 la after 150 jus is 1960.7 m/s, correspond- 
ing to a degree of overdrive of /= 1.02.(Note that this calcu- 
lation used a driver pressure ratio of r p = 70). 

The structure of the detonation wave is similar to that for 
the pq = 0.4 bar case, but the frequencies of oscillation are 
much smaller. For 4> = 1.0, the frequency is 0.34 ± 0.02 MHz, 
and it decreases to 0.29 MHz for (() = 0.8. At (f> = 0.6, the high 
frequency oscillations are much less uniform and range from 
approximately 0.18 MHz to 0.28 MHz. At this equivalence 
ratio ((() = 0.6) the low frequency mode begins to become evi- 
dent. At 4> = 0.4, the only mode present is the low-frequency 
one, except during maximum detonation speed when very 
weak traces of the high frequency mode can be seen. 

The results for p 0 = 1.0 bar are shown in fig. 12. A finer 

grid, having a minimum spacing of Ax Inin =7.324xlO" 5 cm, 
was used in this case, due to the shorter ignition delay times 
at this pressure near stoichiometric conditions. The theoreti- 
cal Chapman- Jouguet detonation speed for a stoichiometric 
hydrogen-air mixture at p 0 = 1.0 bar and T 0 = 298 K is 
1976.7 m/s. The average detonation speed computed in fig. 


12a after 15 ps is 2127.7 m/s, corresponding to a degree of 
overdrive of/= 1.16. Note that this calculation used a driver 
pressure ratio of r p = 150. In general, the higher the pressure, 
the higher the degree of overdrive needed to obtain the 
high-frequency, low amplitude mode of oscillation described 
by the McVey-Toong mechanism. 

The results for p 0 = 1.0 bar show more variation in the 
amplitude of oscillation than the p 0 = 0.4 and 0.2 bar cases, 
even though the degree of overdrive is higher in this case. 
The frequencies of oscillation are also much higher at 
Pq = 1.0 bar. For the 4> = 10 case, the frequency of oscillation 
is 3.2 ± 0.2 MHz. At lower equivalence ratios, the frequen- 
cies are 2.71 MHz for 4> = 0.8, 2.08 MHz for 4> = 0.6, and 
1.28 MHz for 4> = 0.4. Note that due to the high degree of 
overdrive, the high-frequency mode persists at the smallest 
value of the equivalence ratio (4> = 0.4). 

Driver pressure ratio e ffect 

The effect of driver pressure ratio, r p , was investigated 
for hydrogen-air mixtures at p 0 = 0.4 bar. Figure 13 shows 
the detonation speed for four mixtures at a driver pressure 
r p = 150. It can be seen that compared with the results of 
fig. 8, (which were obtained with r p = 100) the oscillations 
are much more uniform, and have a lower amplitude and a 
higher frequency. The high-frequency mode is present even 
for the 4> = 0.4 case. 

At 4>=10, the average detonation speed after 50 ps is 
2087.3 m/s, corresponding to a degree of overdrive/= 1.14. 
The frequency of oscillation for this case is 1.19 ±0.03 
MHz. The frequencies for the other cases are: 1.09 MHz for 
4>=0.8, 0.91 MHz for 4>=0.6, and 0.63 MHz for 4>=0.4. 

Figure 14 shows the results for the same mixtures, but a 
driver pressure ratio of r p = 40. Previous “steady” detonation 

calculations 3,4 (i.e. coarse-grid computations) using this 
driver pressure ratio, have shown that the resulting detona- 
tion speeds agree closely with the Chapman- Jouguet theoret- 
ical speed. That is, for this driver pressure ratio, the degree of 
overdrive should be very close to one (i.e.,/= 1.0). However, 
for the high-resolution computations in the present study, the 
detonation wave at this driver pressure ratio contains very 
complex, irregular oscillations for the equivalence ratio 
range 0.6 < 4> < 1.0. (The 4> = 0.4 case should be run for 
longer times to determine the oscillation pattern). This 
behavior for p 0 = 0.4 bar, and degree of overdrive close to 

one, agrees with the findings of Sussman 18 for a stoichiomet- 
ric hydrogen-air mixture at nearly identical initial pressure 
and temperature. Bourlioux et al . 2J also showed that the 
oscillations obtained with a one-step chemistry model 
became nearly chaotic for low degrees of overdrive. It is not 
clear wether the one-dimensional detonations considered 
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here are inherently unstable at these conditions or wether 
further grid refinement is necessary to obtain accurate solu- 
tions. 

Conclusions 

The one-dimensional Euler equations with finite rate 
chemistry were solved to investigate the structure and stabil- 
ity of detonation waves for hydrogen-air mixtures at various 
equivalence ratios and pressures. The results indicate that the 
solutions to the governing equations do not include a steadily 
propagating detonation wave. Instead, the solutions always 
exhibit an unsteady oscillatory propagation mode. For mix- 
tures close to stoichiometric conditions, the detonation wave 
propagates in a high-frequency low amplitude mode that is 
described by the McVey-Toong mechanism. For very lean 
mixtures, the detonation wave switches to a low-frequency, 
high-amplitude propagation mode similar to the “galloping” 
detonations observed experimentally in near-limit mixtures. 

The transition from the high-frequency mode to the 
low-frequency mode depends on the degree of overdrive. 
High degrees of overdrive tend to stabilize the detonation 
wave and delay the transition to the low-frequency mode for 
lean mixtures. The frequency and amplitude of the oscilla- 
tions are very uniform for high degrees of overdrive. At low 
pressures (p 0 = 0.2 bar), solutions that behave according to 
the McVey-Toong mechanism were possible for degrees of 
overdrive close to one. For higher pressures (p 0 > 0.4 bar), 
these solutions were possible only for degrees of overdrive 
above 1.07. Calculations with degrees of overdrive close to 
one at the higher pressures resulted in complex, irregular 
oscillations. 
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TABLE 1. — H 2 -Air Reaction Mechanism 3 


No. 

Reaction 

<4 

b 

0 

1 

h 2 + o 2 =ho 2 + h 

l.OxlO 14 

0.0 

28197.38 

2 

H + 0 2 = OH + O 

2.6xl0 14 

0.0 

8459.21 

3 

h 2 + o = oh + h 

1.8xl0 10 

1.0 

4481.37 

4 

H 2 + OH = H + H 2 0 

2.2xl0 13 

0.0 

2593.15 

5 

OH + OH = O + H 2 0 

6.3xl0 12 

0.0 

548.84 

6 b 

H + OH + M = H 2 0 + M 

2.2xl0 22 

-2.0 

0.0 

7 b 

h + h + m = h 2 + m 

6.4xl0 17 

-1.0 

0.0 

8 b 

H + O + M = OH + M 

6.0xl0 16 

-0.6 

0.0 

9 b 

H + 0 2 + M = H0 2 + M 

2.1xl0 15 

0.0 

-503.52 

10 

o + o + m = o 2 + m 

6.0xl0 13 

0.0 

906.34 

11 

ho 2 + H = OH + OH 

1.4xl0 14 

0.0 

543.81 

12 

ho 2 + h=h 2 o + o 

l.OxlO 13 

0.0 

543.81 

13 

ho 2 + o = o 2 + OH 

1.5xl0 13 

0.0 

478.35 

14 

ho 2 + oh = h 2 o + o 2 

8.0xl0 12 

0.0 

0.0 

15 

H0 2 4* ho 2 = h 2 o 2 + o 2 

2.0xl0 12 

0.0 

0.0 

16 

h + h 2 o 2 = h 2 + ho 2 

1.4xl0 12 

0.0 

1812.69 

17 

o + h 2 o 2 = oh + ho 2 

1.4xl0 13 

0.0 

3222.56 

18 

OH + H 2 0 2 = H 2 0 + H0 2 

6.1xl0 12 

0.0 

720.04 

19 b 

h 2 o 2 + m = oh + oh + m 

1.2xl0 17 

0.0 

22910.37 


a Forward rate coefficient kj- - /\7^exp(-0/T); units are moles, 
seconds, centimeters, and Kelvins. 


b Third-body efficiencies: 

Reaction 6: H 2 0 = 6.0; Reaction 7: H 2 0 = 6.0; 

Reaction 8: H 2 0 = 5.0; Reaction 9: H 2 0 = 16.0, H 2 = 2.0; 
Reaction 19: H 2 0 = 15.0. 


TABLE 2. — Ignition Delay Time 


(ti g ) for ([>=1.0. 

.T^HSO.OK 

Pi (bar) 

tig (Pis) 

0.5 

50.3 

1.0 

25.1 

10.0 

4.8 

15.0 

11.1 

20.0 

41.0 

25.0 

54.7 

30.0 

60.3 

40.0 

62.3 

50.0 

60.1 

75.0 

52.0 

100.0 

44.9 
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Fig. 1 . Ignition delay for Chapman -Jouguet detonations in H 2 -Air; Po=0.4 bar, T 0 =298 K. 
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Fig. 2. Ignition delay for Chapman -Jouguet detonations in H 2 -Air; p 0 =1.0 bar, T 0 =298 K. 
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Fig. 3. Computational grid at three different times. 



Fig. 4. Schematic of detonation wave initiation. 


(a) 






Fig. 5. Detonation speed as a function of time; H 2 -Air; 
p 0 = 0.4 bar, <|>=1.0; r p = Pdri/Po = 100; D CJ = 1956.2 m/s; 
f= 1.07. 
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Fig. 6. Pressure ratio as a function of distance at different 
times; H 2 -Air; p 0 = 0.4 bar, (f)=1.0; r p = p dri Jp 0 = 100. 




Fig. 7. History of pressure and density (x-t diagram) 
showing the Mcvey-Toong mechanism; ^=150. 




Fig. 8. Detonation speed as a function of time; H 2 -Air; 
Po = 0.4 bar, r p = 100, Ax min =2.441xl0' 4 cm. 
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Time (^s) 



Time (ns) 





Fig. 9. Maximum pressure as a function of time; H 2 -Air; 
p 0 = 0.4 bar, r p = 100, Ax m i n =2.441xl0' 4 cm. 





Fig. 10. Detonation speed as a function of time; H 2 -Air; 
p 0 = 0.4 bar, r p = 100, Ax min =2.441xl0" 4 cm. 
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(a) 



(b) 



(c) 



(d) 



Fig. 1 1 . Detonation speed as a function of time; H 2 -Air; 
p 0 = 0.2 bar, r p = 70, Ax min =2.441xl0‘ 4 cm. 


(a) 



(b) 



(c) 



(d) 



Fig. 12. Detonation speed as a function of time; H 2 -Air; 
p 0 = 1.0 bar, r p = 150, Ax mill =7. 324x1 O' 5 cm. 
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Detonation speed (m/s) Detonation speed (m/s) Detonation speed (m/s) Detonation speed (m/s) 










Fig. 13. Detonation speed as a function of time; H 2 -Air; 
p 0 = 0.4 bar, r p = 150, Ax^^^lxlO" 4 cm. 


Fig. 14. Detonation speed as a function of time; H 2 -Air; 
p 0 = 0.4 bar, r p = 40, Ax^^^lxlO -4 cm. 
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